Quantized spin Hall effect in 3 He-A and other p-wave paired Fermi systems 
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In this paper we propose the quantized spin Hall effect (SHE) in the vortex state of a rotating 
p-wave paired Fermi system in an inhomogeneous magnetic field and in a weak periodic potential. 
It is the three dimensional extension of the spin Hall effect for a 3 He-A superfluid film studied in 
Ref.— . It may also be considered as a generalization of the 3D quantized charge Hall effect of Bloch 
electrons in Refi2 to the spin transport. The A-phase of 3 He or, more generally, the p-wave paired 
phase of a cold Fermi atomic gas, under suitable conditions should be a good candidate to observe 
the SHE, because the system has a conserved spin current (with no spin-orbit couplings). 

I. INTRODUCTION 

Recently, spin transport in condensed matter systems has received wide attention. One of the most striking 
phenomena is the spin Hall effect (SHE) in semiconductors 3 -. In this effect, the spin current is driven by an electric 
field, and flows in a direction perpendicular to it. Time reversal symmetry is respected, the situation is remarkably 
different from the Hall effect of charged currents in a magnetic field. However, the SHE in semiconductors needs 
appreciable spin-orbit couplings, which make the spin current non-conserved. This causes several complications 
in theoretical study, and makes the relationship between the spin current and the experimentally observable spin 
accumulation intricate and obscure. 

In the literature, another class of condensed matter systems, namely supcrfluids, have been suggested in which the 
SHE may be observable without the need of spin-orbit couplings. It has been shown that in the vortex state of two- 
dimensional (2D) d-wave superconductors, the spin current can be generated by an inhomogeneous magnetic field and 
flows perpendicular to the field gradient. Moreover, a similar spin Hall effect was suggested for a rotating superfluid 
3 He film in a non-uniform magnetic fieldi. In this case, time reversal symmetry is broken. If the quasiparticle spectrum 
has a gap and the Fermi level lies within the gap, the spin Hall conductance can be expressed as a topological number 
(the first Chern number) and therefore is quantized, in the units of /iHe/87T with K = c = 1 where //He is the Zeeman 
coupling constant for 3 He atom. This effect is thus the spin analogue of the quantized charge Hall effect of Bloch 
electrons^. The distinctive feature of the SHE in supcrfluids is the absence of spin-orbit couplings. In general, 
spin-orbit couplings exist in superconductors or semiconductors, because of the presence of the internal crystal fields 
or confining fields (in the case of heterostructures). But spin-orbit couplings do not exist in superfluid 3 He, of which 
the constituents are neutral atoms. (Recall that spin-orbit couplings exist only for charged particles). Superfluid 3 He 
is in the spin-triplet p-wave stated and the spin rotational symmetry is broken spontaneously. However, one of 3 He 
superfluid phases, i.e. the 3 He-A phase, has a residual rotational symmetry around the z-axis in the spin space. In 
this phase the system has a conserved spin currenfeii Therefore, 3 He-A film should be a good candidate to observe 
the QSHE without worrying about complications due to non-conservation of the spin current. This motivates us to 
study this system in more general three dimensional cases. 

In this paper, we propose the quantized SHE in the vortex state of a rotating p-wave paired Fermi system in a non- 
uniform magnetic field and in a periodic potential in three dimensions. Our proposal is a generalization of the SHE 
for a thin film of 3 Hc-A superfluid studied in refi On the other hand, it is the spin analogue of the three-dimensional 
quantum charge Hall effect (3DQHE) for Bloch electrons suggested in refi^. Our discussions are valid also for the 
vortex state of a cold Fermi atomic gas that is in the p-wave pairing regime^. 

II. VORTEX STATE IN THE 3 HE-A PHASE 

We use the same notations as in refi The gap matrix for a 3 Hc superfluid is defined as A Q ^(x, y) = 
J2fS ^a/37<5( x J y)(' t P-y(' x -)' l Ps{'y)) , where ip a (x) is the Fermion field for the 3 He atom with spin a =|, |,(...) the thermal ex- 
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TABLE I: The correspondence between a superconductor in an applied magnetic field and a rotating superfluid. 
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pectation value, and V a /3 7 s(x, y) the pairing interaction. Under the U(l) gauge transformation i/j a (x) — + e _ie ( x ^ a (x), 
the gap matrix is transformed as A Q/ 3(x, y) — > A a( g(x, y)e _ie '- x ^ _ie ^ y ^. The matrix can be written in a form 
A Q/ g(x, y) = $ (R(?) / (2ot)3 (k)e lk '( x ~ y ) , where $ (Rg) denotes the gap amplitude depending on the center of 
mass of a Cooper pair R G = (x + y)/2 and ^ Q| g(k) the order parameter. For the spin triplet pairing state, the order 
parameter is written by using d(k)-vector in the spin spaced as \fr a fl(k) = id(k) • (<J<J2) a p, where a = (ai, 02, 03) is 
the Pauli matrix. 

In the 3 He-A phase, the d-vector is given by 

d(k) = e z (k x + ik y ). (1) 

Clearly the order parameter is invariant under rotations around the z-axis in the spin space. As we mentioned, 
the presence of the spin rotational symmetry is crucial for defining a conserved spin current. In the particle- hole 
representation, the quasiparticles arc described by a two-component field 

The Hamiltonian for a rotating superfluid with angular velocity Q in the rotating frame readsi 

H = J d 3 xdW(xMp,x,y)*(y) (2) 



with 



/ e(p-mR)5 3 (x-y) A£ v (x, y) e -^( x )+^(y)) \ 
^ P ,x, yj ^ -Af (x,y)e5^W+^y» -e(p + mR)5 3 (x - y), ) 

e(pimR) = — — (p ± toR) 2 — up, 
2m 

p = -iV, 

r = n x x, 

where A™ v (x, y) is the gauge invariant part of the spin-traced-out gap function for A-phasc Aa (x, y) = 
(l/2)Tr[A(x, y)cri] and (y(x) + ip(y))/2 the gauge dependent phase. Namely, we have a relation Aa(x, y) = 
A A lv (x, y) exp[— i(<p(x) + tp(y))/2]. Moreover, since A^ v (x, y) depends only on the relative coordinate x y, for 
p-wave pairing we have A A lv *(y,x) = — A A lv *(x, y). From the above equations the correspondence between the case 
of a superconductor in an applied magnetic field and that of a rotating superfluid is clear, as shown in Table |TJ From 
the correspondence, one naturally expects that a vortex state emerges when Q is larger than a certain critical value 
f2 c i. The vortex lines are parallel to fi, and the vortices form a lattice on the plane perpendicular to n = Jl/O. The 
center x l of the i-th vortex can be expressed in terms of the primitive vectors a and b of the lattice as 

x' = Z^a + njb + zn, (ij, rii integers; z real), 

where the phase </?(x) of the gap function satisfies 

VxV</?(x) = 27rn^](5 3 (x-x l ). (3) 

i 

Since a Cooper pair has a "charge" 2m, a vortex has a "flux" ir/m, and one obtains 

„ TT 

to a x b 
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This relation is consistent with the correspondence between the fi-flux in the unit cell of the vortex lattice and the 
flux quanta ir/e in superconductors which is a half of the unit flux 2n/e. The difference is due to the fact that the 
charge of a Cooper pair in superconductors is 2e. 

Now we proceed to examine the lattice symmetry of this vortex state. We will see that the system has a " magnetic" 
translational symmetry with respect to two cells of the vortex lattice, since a vortex has a "flux" n/m, which is a half 
of the "unit flux" 2ir /m. In an analogy to the system of charge particles on a lattice, we consider the "magnetic" 
translation operator- 
Tar = exp[i(5r • (p + t 3 toR)], 

where t 3 is the third component of the Pauli matrices in particle-hole space. By using the gauge degrees of freedom, 
one can introduce the following constraints without conflict with Eq. |5]h 

cp(x + a) = </?(x) — ma • R, 
y(x + b) = ip(x.) — mb • R. 

Then we see that the system has the translational symmetry generated by T a and T 2 b 

[H,T a ] = [H, T 2b ] = 0, [T a , T 2 t>] = 0. 

Let us suppose that there is a periodic potential, e.g. as what happens in an optical lattice 

U{x) = [/ cosK-x. (4) 

Then the Hamiltonian density becomes 

H(p, x, y) Hu(p, x, y) = H(p, x, y) + C/(x)t 3 <5 3 (x - y). (5) 

Suppose that K is not in the a6-plane, i.e. K-n^0 and consider a cell spanned by the following vectors 

a' = a + a n n, 
2b' = 2(b + /3„n), 

n, (6) 



Kn 

where 

f 



(2tt-K- a), 



K n 

Pn = K^^- K ' b )- 
We see that our system has a three dimensional periodicity represented by the commutation relations 

[H,T aL ,] = [H,T 2b ,] = [H,t c ,}=0, (7) 
where t c * = exp[ic' • p] is the usual translational operator. We also note that 

[T a ,,T 2b ,] = [T 2V ,tcf] - [tcfM = 0. (8) 
The Bogoliubov-de Gennes equation for our system reads 



J d 3 yHu(p,x,y)<S>E(y) = E<P E { X ) 

-Vg(x) 



Mx) = I ) , 



where ([/^(x), Ve(x)) is related to the fermion operators as 
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TABLE II: The mapping between physical quantities in the Hamiltonian for Bloch electrons in QH systero^^ and that for 
superfluid 3 He in the particle-hole representation 



Bloch electrons 


superfluid 3 He 
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1/2 


Pe 


Pi 


Je 


y a 
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VB 3 



and r y Ea and ^Eu are the creation and annihilation operator, respectively, of the Bogoliubov quasiparticles. In 
accordance with Eqs. ([7]) and ([8]), the solution of Eq. (|9|) is of the Bloch form 

$nk(x) = e ik ' x w, lk (x), 

E = -E ra k, 

where n is the band index and k is the crystal momentum defined in the 'magnetic' Brillouin Zone 

—ix/a! < k x < 7r/o', 
-it /2b' <k y < 7r/2&', 
— 7r/c' < k z < n/c', 

and u n k(x) is commensurate with the 'magnetic' unit cell. Here we assume the existence of a gap in the quasiparticle 
spectrum around zero energy. 



III. 3D SPIN HALL CONDUCTANCE AND CHERN NUMBER 

We define a spin current from the spin conservation law 

p s (x) + V-j* = 0, (10) 

where 

is the density of the z-components of the spins. The conservation is assured by the spin rotational symmetry around 
the z-axis in the 3 He-A phase Eq.(fT]). Note that, in the particle-hole representation, the spin density is written by 
the sum of $^ k (x)$ n k(x) which has the same form as the charge density of a state (n, k) in the usual representation. 
From Eqs. (fit))) and (jTTJ) , the explicit form of the spin current for the z component is given by 

which is also of the same form as the charge current in the usual representation. 

To drive the spin current, we apply an inhomogeneous magnetic field in the z-direction, B = B z e zi with a constant 
gradient, i.e. VB 2 = const. Then to the Hamiltonian (|5|) the Zccman coupling term TLjj should be added, which in 
the particle-hole representation reads 

H I = ^x-VB z . (12) 

This has the same form as the interaction between the charge density and a constant electric field in the usual 
representation. Then, one can see a mapping between physical quantities in the Hamiltionian for Bloch electrons in 
the quantum Hall (QH) system™^ and that for superfluid 3 He in the particle-hole representation. See Table HT1 

The mapping tells us that the expectation value of the spin current obtained by the Kubo formula^ has the same 
form as that for the charged Hall current of 3D Bloch electrons 2 ; namely, 



<Jti>=< j V j B z , 
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with a conductivity 

f d 3 k < Unk\Jsi\u m k>< U m k\J"j\Unk > —{u <-> m) 

LlBZ (Ernk - E n ^Y 



n < . m 



= ^E/ 7|t^[V k xA n (k)] fe) (13) 

where Vk = (9/<9k and A(k) is a gauge connection in k-spacc 

A„(k) =< u nk |V k |u, ik >, 

which was first introduced in this context by one of the authors^. 

Following Refi^, we can relate the spin Hall conductance (|13[) to the topological Chern number— £ as follows. The 
crystal momentum is parametrized as 

k = f a ,G a ,+f b ,G b ,/2 + f c ,G c , (14) 

with < fa',.h'Jc' < 1, and 

f At' Gc ' ' C ' t j2i„ _ 2lT f 1 jr 



d A k = / df' c — — — / d z k = — dfcf / d^k, 

MBZ JO c> JT 2 (/ C ,) c ' JO JT 2 {f c ,) 

where T 2 (/ c <) is a two dimensional torus formed by the k's with fixed f c >. Similarly, 

/ d 3 k = ^ [ df a , [ d 2 k=?- f dfy [ d 2 k. 

JMBZ a JO JT 2 (f a ,) Jo JT 2 (f b ,) 

One may introduce a vector D that satisfies 

&ij = CijkDk- (15) 
We write D in terms of the primitive vectors of the reciprocal lattice as 

D = aG ,+/3G 6 //2 + 7 G c ,, (16) 

and obtain 



7 27T- 



lc'-D = -^ f dfc ,^2 I £4 [Vu x A„(k)] • ^ 

2tt 16tt 2 7 ^ o J T 2 {fcl) 2m \c'\ 

Therefore, the integral 

W p: [V k X A„(k)] • * 



(17) 



is the first Chern numbes^ on the torus in momentum space. For the discussions of the relation between Berry's 
phased and Eq. (jTTJ) , see Refsi^. Because of its topological nature, the integral is an integer which depends on f c i 
when the Fermi level for the quasi-particles is in a gap between magnetic Bloch bands. Then, 

7 = -T^', (18) 
Ion* 



with N c i an integer. We can also show that 



13 = -T^ Nb '' (19) 



where N a > and arc integers. So finally we have 

MHe 
16V 2 

where G is a lattice vector reciprocal to the original lattice formed by a',2b',c'. Thus, the spin current is shown to 
be 

I6n~ 2 

This result is the spin analogy in the present system of the charged Hall current in 3DQHE 2 



D = -tF2 G - (2°) 



< J >=SGxVB, (21) 
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IV. SUMMARY AND DISCUSSIONS 



In summary, superfluid 3 He-A is a good candidate system to investigate the SHE, since this system has no spin-orbit 
couplings, and has a spin rotational symmetry which ensures the spin current conservation. We have presented the 3D 
extension of the quantized spin Hall effect in the vortex state of a rotating 3 He-A superfluid filmi in an inhomogeneous 
magnetic held. We have demonstrated the one to one correspondence between the SHE in 3 He-A in the case at hand 
and 3DQHE of Bloch electrons previously studied in refi£, and thus have predicted the existence of the 3d QSHE in 
this system, when suitable conditions are satisfied. 

Here we summarize in more details the conditions that are needed for the predicted QSHE. As we have mentioned, 
the 3 He-A superfluid has to be in the vortex phase. The vortices normally form a regular triangular lattice with the 
primitive vectors a = (a, 0, 0) and b = (a/2, V3a/2, 0) in the plane perpendicular to the angular velocity ft = (0, 0, fl). 
Then we apply an one-dimensional periodic potential, as what happens in an optical lattice, along an axis which is 
not in the plane of the vortex lattice (the plane formed by vectors a and b. In the simplest case, this axis can coincide 
with the axis of rotation, i.e. that of fi. This periodic potential should be weak, not to spoil the superfluidity and the 
vortex phase. Then to drive the spin current we need an inhomogeneous, static magnetic field in a fixed direction, 
coupled to the nuclear magnetic moment of the 3 Hc nuclei, with constant gradient in its magnitude. When these 
conditions are satisfied, we have shown an exact analogy between the spin current of the quasi-particles in the present 
system (in the particle-hole representation) and the charge current of the Bloch electrons in a three dimensional 
periodic potential. Based on the knowledge of latter— with a connection to the Diophantine equations (see ref.— and 
the third paper in Ref £) , a three dimensional quantized spin Hall effect with non-zero Hall conductances is generically 
expected, when the Fermi level of the Bogoliubov quasi-particlcs just completely fill a certain number of Bloch bands, 
opening a gap in the spectrum of quasi-particles. 

Since the vortex lattice is created by rotation, the vortex structure near boundary could be rather complicated. How 
spin accumulation would be affected by conditions at the boundary is an important issue for experimental detection. 
Theoretically this is a separate issue from the bulk effect we have addressed in this paper, and is certainly worth 
further study. 

Finally we emphasize that in the above we have mentioned the 3 Hc-A phase as a familiar example of a p-wave paired 
Fermi system. Actually the above discussions are valid also for cold Fermi atomic gases that allow p-wavc pairing in 
a certain parameter regime. This remark opens the door for more candidate systems that may experimentally realize 
the predicted QSHE. 
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